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Abstract 

We derive the classical continuum limit of the operator tr s a tJiV D ov {x, x) with 
D ov being the overlap Dirac operator and show that it corresponds to the gauge 
field strength tensor F^ u (x). 
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1 Introduction 



Lattice gauge operators are usually constructed explicitly from link variables. For 
example, the Wilson gauge action uses the product of gauge links at the boundary of 
the square plaquette. Similarly, the gauge field strength tensor can be defined through 
suitable combinations of such elementary plaquettes. To improve scaling behavior of 
the action and other gauge operators, rectangular and more sophisticated loops have 
been included [U [21 El H] . Furthermore, it was shown that operators with smeared 
gauge links, being less ultralocal, are effective in filtering out ultraviolet fluctuations 
and improving the efficiency of numerical calculations [3, El [71 18] . 

A different approach is realized in constructing topological charge densities via 
Ginsparg- Wilson Dirac operators It uses matrix elements of the Dirac operator as 
a starting point. By virtue of the fact that this type of Dirac operators are inevitably 
non-ultralocal [10], they entail a sum of all gauge loops and are thus automatically 
smeared. While technically involved, this approach has the advantage that Ginsparg- 
Wilson operators incorporate an exact lattice chiral symmetry [IT] . This leads to an 
index theorem on the lattice [9]. 

The overlap Dirac operator [12J offers a concrete example which satisfies Ginsparg- 
Wilson relation [13J. Its compact form makes it amenable to analytic as well as nu- 
merical calculations. The associated topological charge density was shown to have 
the correct classical continuum limit via weak-coupling expansion [TJ] and by direct 
calculations [T51 [TBI E]- It is with this topological density operator that the sub- 
dimensional long-range structure has been discovered [TS1 [TjJl [20] and confirmed [2JJ 
in QCD, as well as in 2-D CP(N-l) models [22]. It is also with this operator that the 
required negativity of the topological density correlator is borne out clearly with only 
a handful of configurations [23J. Whereas, detecting this negativity to the same preci- 
sion using conventional operators, such as those used in the glueball calculation [24], 
would require much larger statistics. Based on these observations, it was suggested 
by one of the authors that the condition of chiral symmetry plays a relevant role in 
efficient suppression of the ultraviolet noise [25], and that all gauge operators can 
be constructed from the chirally symmetric Dirac operator [26J. This way, one can 
also have a formulation of lattice QCD where the gauge action, the 6 term and the 
fermion action are all expressed in terms of the lattice Dirac operator [25] . 

In the present work, we will concentrate on the gauge field strength tensor. It was 
suggested that the classical limit of the tensor component of the overlap operator, 
tr s o^y D ov (x, x), is proportional to the gauge field strength tensor F^ u (x) [261 I2~T] . 
We shall explicitly calculate it and show that this is indeed the case. In view of the 
fact the chirally symmetric Dirac operators are non-ultralocal [10], the derivation is 
somewhat non-trivial. Since the overlap Dirac operator is expected to be local on 
gauge configurations of interest [281 [29], our result implies that the quantum operator 
constructed this way represents a valid definition of gauge field strength tensor in 
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lattice QCD. This could provide a practical tool for evaluating gluonic observables. 
A preliminary version of the calculation was given in Ref. [30] . 

2 Formulation 

Our goal in this paper is to show that if we discretize a classical gauge field A^(x) on 
a hypercubic lattice with (classical) lattice spacing a, and consider the overlap Dirac 
operator D ov on such background, then 

tr s D m (x, x) oc a 2 F^ u (x) for a — >0 (1) 

where "tr s " denotes the spinor trace. Q By classical SU(3) gauge configuration we 
mean any configuration of the gauge field that is smooth (differentiable arbitrarily 
many times) almost everywhere. In the above equation we have implicitly assumed 
that x is not a singular point of A^(x). It was also implicitly assumed that point x 
is a lattice grid point of a superimposed lattice for arbitrary lattice spacing a. This 
is technically most easily (and without loss of generality) achieved if we focus on the 
point x = with the origin of the lattice coordinate system aligned with the Cartesian 
coordinate system of R 4 . Also, to avoid an extensive discussion of technical issues 
associated with the transcription of the field with singularities onto the hypercubic 
lattice, we will implicitly assume in the following that the field is smooth everywhere. 
This is sensible since due to the locality of the overlap operator, the result is expected 
to be valid for an arbitrary non-singular space-time point x. For a discussion relevant 
to this point, the reader is referred to Ref. [15] . 

In the convention that we will use, the continuum gauge potential A^(x) is the 
vector field of traceless Hermitian matrices and the corresponding field-strength 
tensor is 

Ffiv(x) = d^A u (x) - d u A^(x) + i[A^x), A u (x)}. (2) 
With the covariant derivative defined as 

D^ix) = (d^ + iA^x))^), (3) 

one has 

[D„,D v \^{x) = iF^{x)(j){x). (4) 

1 Note that we use the convention that the real-valued arguments of lattice quantities (such as 
D(x,x) on LHS of Eq. ([1])) are given in parenthesis, while the integer-valued lattice coordinates are 
written as subscripts (such as D n ^ n ). 

2 Note that this differs from conventions of Ref. [55], where anti-Hermitian gauge potentials were 
used instead. The equations below can be obtained from equations of Ref. [26) via substitutions 
A fl (x) — > iA^(x), F liJ/ (x) — > iF^ v {x). The value of constant c T in Eq.© is the same in both 
conventions. 
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The transcription of A^x) to the hypercubic lattice with integer coordinates 
n = (ni, n 2 , n 3 , n±) is accomplished in a standard manner. If a is the classical lat- 
tice spacing, we associate the lattice site n with the space-time point x = an, and 



the lattice link variable U ntll is defined as 



U n A a ) = exp(mA M (<m)) (5) 



The overlap Dirac operator D ov is given by p2] 

D" = ph+X-j===y t X = y-R-p + 4r (6) 

where — p, p G (0,2r), is the negative mass parameter and 

D, = \{U,S, - SjflJ] R= T - Y\U»S„ + fijflj] (7) 

with 

(<S'/i)m,n = 3m,n—fi {U ^i)m,n — U m ^ §m,n (8) 

We shall take the Euclidean 7-matrices to be Hermitian, i.e. 7^ = 7^ and {7^, 7^} = 
2<W- 

With the above defining relations, we shall proceed to show the following in an 
explicit calculation: 

If Afj,(x) is a smooth SU(3) gauge potential on i? 4 , and U(a) is the transcription of 
this field to the hypercubic lattice with classical lattice spacing a, then 

tT s a^DH(u{a)) = c T a 2 F^(0) + 0(a 3 ), (9) 

where D™ Q is the matrix element of the overlap operator at (to, n) = (0,0). The 
non-zero constant c T = c T (p) is independent of A^x), and a^ u = ^[7^,7^]- 



3 Calculation 

To proceed with the calculation, we shall assume that tr s a^ v Dq v has a Taylor ex- 
pansion in a and we will compute the leading contributions. 

To evaluate the diagonal element D™ n we introduce the momentum variable in 
the following way [T5], [T7] 



D 



PIT fj^h 

V D° V 5 = V D° V / -ikjn-m) QqN 

m m J ~ 7T l Z7r J 

' —ikn \ "> 7~)cw ikm (^^\ 



n (2vr) 4e 



m 



This will allow us to evaluate the inverse square root in Eq. ([6]). Next, we define the 
diagonal matrices K(k) as 

(K(k)) n , m = e ikn 5 n , m . (12) 

These matrices are unitary: K^(k) = K~ 1 (k) = K(—k). If we now introduce the 
vector 1 such that l n = 1, a vector with all entries set to 1, we can rewrite the above 
expression as 

To calculate K~ 1 (k)D ov K(k), we assume that we can express as a P ower 

series in X^X. Then 

K-\k)D™K{k) = p(l- y^=) , (14) 



where 
and 



Y = -X = -K-\k)XK(k) = M + R-ty-ift, (15) 



i(e^(r7^-l)-e-^(fij^-l)), 
R = r :;E{^(U,S»-l) + e- ik »(SlUl-l)] 

M = p + r5>A-l), (16) 



where s M = sin/c^, = cos/c M . 



Eqs. (jlDjl . (TBI), and (ITD. lead to 

r* d 4 fc 1 



tr ,^ flsr „ = _,£_ tr ,^ (k^i) . (17) 



3.1 Computational strategy 



As we mentioned before, we assume that tr s a^ v D™ has a Taylor expansion in a. 
We will compute the leading contributions by taking derivatives with respect to a 
and then evaluating the limit a — > 0. We assume that all the matrices and matrix 
products are well defined and that we can take derivatives in the usual fashion. The 
non-trivial part of the calculation is taking the derivative of -7=yp with respect to a. 
For this we will write 

1 r°° , l 
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under the condition that Y^Y does not have zero eigenvalues, which is satisfied for 
the case of classical gauge fields sufficiently close to the continuum limit [28J . Noting 
that (M- 1 )' = -M~ l M'M- 1 for a matrix M, we have 

1 . ijjW 1 



d 2 1 2 /•«> 1 /d.. t .\ 1 / d , \ 1 



rfd r— —Y^Y — — — F f F 



da 2 v / yty ttJ-oo cr 2 + W \^da y <j 2 + Yty \^da / a 2 + yty 

1 ™ i / H 2 \ 1 

' ^ = ^7^^ „ =— ; • (19) 

ttj-oo a 2 + yty\ v da 2 ;a 2 + yty v ; 

We see that the problem is reduced to computing derivatives of Y with respect to 
a. The only matrices that depend on a are the link matrices U^. Since we are only 
interested in the limit a — > and since we will only calculate the contributions up to 
order a 2 , we only need the following limits 

hmU, = 1, 

limJ-U, = L4 M (0)1, 
a-»o da 

lim-^E^ = -A lx (0) 2 l + 2iNVA lx (0). (20) 

where 

(JWA M (0)) m , n = E^^(°)^,n, (21) 

with n a being the component of the 4-vector n. We see that in the limit a — > both 
f/ M and j-L^j reduce to an identity matrix in the space-time coordinates, but they are 
not necessarily diagonal in color space. However, the second derivative NVA^O) has 
a term that is different: this matrix is still diagonal in the space-time index (since it 
is the derivative of a diagonal matrix) but the diagonal elements are not equal. 

To compute (F^=l)o we will need to justify several relations. 

• Relation 1: 



where 



YjY l = zl, (22) 



Y = lim Y = M + R + % + ifi, (23) 

a—>0 

= limD, = i(e^(^-l)-e-^(^-; 
Ro = limi?=^(e^(^-l) + e-^(St-l)), (24) 

a M 
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and 

z = E4 + M2 > (25) 



ft 



is a number. 

To show Eq. (!22j) . we write 

FoVo = ((M + i? ) + (P + ifi)) (( M + ^o) " (Pb + ^ 

= (M + i^o) 2 - (j& + ifi) 2 + [^ + ifi,M + ^o]. (26) 
The commutator is zero since [S^, S u ] = and thus we have 

ySYq = M 2 + R 2 + 2MRq + fi 2 - - *{fi, r> }. (27) 

It is easy to see that RqI = and JZ) 1 = since the "derivative" like term, — 1, 
vanishes when acting on 1. This "derivative" term comes from the continuum limit 
of {U^S^ — 1) which becomes (S^ — 1) as a — > 0. We thus have 

yJy 1 = (M 2 + 8 2 )1 = zl. (28) 

• Relation 2: 

lim - — r- 1 = —r— l = -z - — 1, (29) 



This can be straight-forwardly shown if we expand 2+y ^y as a P ower ser i es m 



and apply Relation 1 in Eq. ff22l) successively. 
• Relation 3: 



—L^JWVO)! = ^^(0)1 - ^-^AiVV^IO)!, (30) 

where 

A = 2Mi? - z{B, %} = 2MR - 2i £ ^(£ M ) . (31) 

To compute the second order derivative in Eq. (|T9|) we will need this relation when 
the matrix - 1 . — acts on non-constant vectors of the form N\7 A U (0)1 in Eq. (1201) . 

<H+Y 'Yb 

To show this we write YqYq as 

FoVo = z + Rl - I/) 2 , + 2Mi? - *{fi, (32) 

V v ' V v ' 

A 2 A 

where A2 and A are unrelated. Using the fact that 

(S v - l)jVVA„(0)l = ((JV + i>)VA M (0)-JVVA M (0))l 

= ^^(0)1 = 3^(0)1, (33) 
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which is a constant vector, we can easily see that A2iVV A M (0)1 = and AAiVVv4 M (0)l 
0. This is because these terms include double "derivatives" like (S^ — 1)(S U — 1) and 
since the first "derivative" produces a constant vector, the second "derivative" act- 
ing on it makes it vanish. Now that we have YqYqN'VA^I — {z + A)NVA^1 and 
AANVA^(0)1 = 0, we can prove by induction that 

(Y^Y ) k NV A^l = z k NVA fl l + kz^ANVA^l. (34) 

From a series expansion P(YqY ) = - 1 1 — and Eq. (I3~4"j) . we see that 

<y Vo 

P(Y^Y )NVA^1 = P(z)NVA li l + P'(z)ANVA fl l. (35) 

Since P(z) = and P'(z) = — , a i\ z ^ , Eqs. ( |34l) and (1351) lead to Relation 3 in 
Eq. dSD]). 

As we mentioned before, we assume that our function admits a Taylor expansion. 
In this case, we can write 

tr s o~nv \ Y . =1 ) = limtr s a uu \ Y , =1 ) 



+ a lim tr s cr^ — ( Y—^==l 



+ ^Si'^-^^TFTK 1 )/ '" 3 '' (36) 



We will now proceed to carry out our calculation order by order. 



3.2 Calculation details 
3.2.1 Order 

We need to compute 

lim tr s [Y-==l) =\ \ da tr s a \Y — 1 1) , (37) 

where the relevant term is tr s a^ v ^o ^+yty 1 = ^r^tr s Y 1. Now 

F 1 = (M + Ro-i8-y> )l = (M-ifi)l, (38) 

due to the fact i? l = I/) 1 = 0. Since Y 1 has only scalar and vector spinor com- 
ponents and no tensor component, it leads to t^o^lo! — and the zeroth order 
contribution is thus zero, i.e: 



s 



3.2.2 Order 1 



We need to compute 

d ( 1 \ 1 f°° f 1 \ 
limtr s a uu — Y - 1 = — / da tr s cr WJy F n ' 1 — 1 

a^o ^daV vW /o 7r -'-°° V <r 2 + y t >o /o 

--/ dtrtr. o, v y (ytyy 1 . (40) 

The first term contribution vanishes since 

Yo 1 = -T—Y'l = -^—{R' - (41) 

and then we are left again with only scalar and vector spinor components. Thus 

To evaluate the second term we need to compute (Y*Y)' l = (2MR' -2iJ2 a s a D' a0 )l 
using the following identities 

K = ^£*AX0)(e^-e-^t), (42) 

D',o = liA,(0)(e^S, + e-^Sl). (43) 

We have then R' l = — r J2 a s aA*(0)l and D' aQ l = ic a A a (0)l and 

(YiY)' l = 2 ^ s a A a (0) (-Mr + c a ) 1, (44) 

a constant vector. We have then 

1 



V + FoVo^V + yoVo 1 



-±—2 £ *«4,(0) (-Mr + c a ) Vo-, x . 



^ 2 + YJY 



' 1 2^s a A a (0)(-Mr + c a )y l 

2 



a 2 + 2 



= (-2^—) 2 ]T s a A a (0) (-Mr + c a ) (M - ifi)l. (45) 

We see that we only have scalar and vector spinor components and this term vanishes 
too after taking the spinor trace. Thus the first order contribution vanishes. 

H^i^vfw 1 ) =°- (46) 
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3.2.3 Order 2 



The derivation of the second order contribution is somewhat involved, but employs the 
same steps as above. The details of the calculations are presented in Appendix A for 
the perusal of interested readers. The main result is that the second order contribution 
is not zero and is given by 

- d A k 
- — — hm tr s a 

n (27r) 4 a-o ' _ v v , , , / u 

/■* d 4 k A{Mc^c u + rs 2 c v + rslc^ 
= (2^ ^ (0) - (47) 

Together with the (null) results from the zeroth and first orders in Eqs. ([39]) and 
the final result is 




tr.^SJ = „^(0),£ § f {M ^ + ^ + 0(A (48) 

Comparing with Eq. flU]), we find 

° {P) = P L J2k? ^ • (49) 

With r = 1 and p = 1.368 (which corresponds to k — 0.19 in the Wilson Dirac 
operator), we find c T = 0.11157.1 



4 Conclusions 

We have shown in an explicit calculation that, for the overlap Dirac operator, the 
classical continuum limit of tr s D ov (x, x) is proportional to the gauge field strength 
tensor F tiv (x). 

Based on the experience of studying the QCD vacuum structure with the topolog- 
ical charge density defined from the overlap operator, it is found that one can obtain 
clear signals with only a handful of gauge configurations [18] HH [20j [21] [22J, [23] . This 

3 We should remark that and R in Eq. ([7|) can be written as D^x) — ^[U fl (x)e adf ' — 
e- ad »Ul(x)] and R(x) = § T,^^)^ + e~ ad "Ul(x)], such as defined in Ref. 02]. Upon taking- 
derivatives with respect to a in e ±adl * and in Eq. (|36|) . the results in Eqs. (|39|) . (J46J) , (|47|) . and 
d48| were obtained [30]. However, due to the fact that in Ref. [30] is defined with an opposite 
sign from the one used here which is <r^„ = ^h^Tv], the result in Ref. |30j is negative of that in 
Eq. 0. 
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is presumably due to the non-ultralocal nature of the overlap operator which serves 
as an efficient filter of the ultraviolet fluctuations [251 126] . It is worthwhile then to 
study whether other operators defined in a similar fashion share this property. For 
example, it would be interesting to see if the calculation of glueball masses, glue mo- 
mentum and angular momentum in the nucleon, etc. can benefit from employing the 
overlap-based definition of the field strength tensor 



which is properly normalized. We should point out that the value of the constant 
c T = c T (p) depends on the mass parameter p used to define the overlap operator. We 
will study this p-dependence in detail elsewhere [3Tj . 

Finally, we wish to mention that, for the purposes of studying QCD vacuum 
structure it is useful to be able to expand gauge observables in low-lying Dirac eigen- 
modes. Indeed, such expansions in the case of overlap-based topological density 
proved to be useful in studying the low-energy behavior of the topological vacuum 
structure [321 EES]- Thus, one rationale for defining all gauge operators in terms of 
Dirac kernels is the fact that it allows such an expansion for an arbitrary opera- 
tor [26]. We note that for the purpose of eigenmode expansion, the expression for 
the field strength tensor in terms of the squared lattice Dirac operators was also 
considered in Ref. [331- 




(50) 
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Appendix A 



We need to compute 



lim tr, a 



s w nv 



Y 



da 2 V v/yty J Q 



datx s a^ Yq 



7T J-oo 
2 roo 



1 

7T J-oo 
1 

a 2 + Y]Y 



datiscr^ Y 



° 2 + YjY J 



° 2 + YjY J 



2 f°° { 
+ - / dat^a^ Y 

7T J-oo \ 



Yn 



° 2 + yJy 
1 



(Y j Y)' 



1 



7T J-oo 



o 2 + YjYo 



i \ 



a 2 + FjFt 



o/o 



The first term in Eq. (1511) turns out to be zero, since the expression 



Y" 

1 n 



Y"l 



a 2 + Y ( jY (T 2 + z a 2 + z 



(51) 



(52) 



only has scalar and vector spinor components and, as a result, it vanishes upon taking 
the the spinor trace in Eq. (|5~T|) . The second term in Eq. (|5T|) also vanishes, because 



Y' 



1 



a 2 + YjY 



(Y^Y)' 



a 2 + Y$Y \a 2 + 



2j2s a A a (0)(-Mr + c a )Y ( ;i (53) 



also has only scalar and vector components. The third term in Eq. fl5Tl) vanishes for 
the same reason, since 



Y 



1 



o 2 + YJY 



i — (Y^Y)' l ——l=(— - — ) Uy s a AJ0)(-Mr + c a )\ Y 1 

/V ,0 \ a 2 + Y ( jY W 2 W V « a> ) 



(54) 

has only spinor and vector components. 

The only non-zero contributions come from the fourth term in Eq. ( }5"TT) . To 
compute its contribution, we need to evaluate 

(ytyyji = ( y ty " + Y »% + 2Y*Y>)1. (55) 
We will compute each term separately. For the first term, we have 

^o" = K - = ~\ E (e*"(7, - r)U; Q S, - e~^( lfl + r)SXo) , (56) 

where f7'' = — A^O) 2 + 2iNVA fl (0). In this case, YqI will have a constant vector 
part 

Wl)e = E (^(0) 2 (i^ - rc M ) + ie-^( 7/ , + r)^(0)) 1, (57) 
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and a non-constant part 

(F Q "1)„ C = -2i£ ArV^(0)(c MT/i - ira^l. (58) 

In considering the first term in Eq. (|55|) . we note that the derivative terms in Y , 

i.e. i?J + IZ> , vanish when acting on constant vectors as shown in Eq. ( |38l) . Further- 
more, the derivatives acting on the non-constant part produce constant vectors as in 
Eq. ( 1331) . As a result, we get a constant term and a non-constant term 

{Y$Y£l) c = (M + ^)^[^(0) 2 ^ 7/ ,-rc M ) + z e -^(7 M + r)a^]l 

- 2i^(c /i 7 M + irs jU )(c 1 ,7 I , - irs,,) 9^^(0)1, 
Mi 1 ' 

(YJn"!)™ = -^M + iBl^W^OKc^-ir^l. (59) 



From Eq. (Job]) , we have 

^ = + $ = £ E (^(7, + - e^( 7 , - r)^^ 1 ) . (60) 



2 , 



Consequently, the second term, Y^YqI, gives 

(Y "V 1) C = - £ [A m (0) 2 (^ 7m + rc M ) + ie- ifc "( 7/ , - r)^(0)] (M - zfl)l, 
(yJ'Vol)^ = 2iE7VV^(0)(c / ,7 /i + ir S/1 )(M-zfi)l. (61) 

The last term to evaluate is 2Y / ^Y '1. We note that since this term only involves 
first derivatives it will only produce a constant vector. Using 

n = ~\ E ^m(O) (e*"( 7/1 - r)£„ + e*>{ llk + r)Sj) , (62) 

we get 

2Y ,t Y 'l =2X)^(0)^(0)(c^ + irs^c^ - irs u )l. (63) 



Mi 1 ' 



Putting all the contributions from Eqs. fl59|) . flBTT) . and fl63l) together, we obtain 
((ytyyji) c = ^[^(0) 2 (-2Mr^-2 S J) + ze-^a^(0)(2Mr + 2z S/ ,)]l 

+ 2 £(^(0)4,(0) - id^A„(0))(c^ + irsjiw - irs v )l, (64) 
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and 

((yty)' 'l) nc = 4^^(0)^(0, - Mr)l. (65) 

We now return to evaluating the last term of the second order contribution in 
Eq. (ELD 



I I I ,. I 

o 

1 \* 1 1 



Yo — j—(Y^Y)" — 1 = -^—Y L_(yty)//i 

V + ltfV ,0 a* + YjY a 2 + z V + *tf V )0 



^ + gJ Y ((Y^1) C + —^—^((yty)^ ( 66 ) 
From Eq. fl30|) . we find for the non-constant term contribution 

Y 2 * ((Y^l) nc = -^—Y ((Y^l) nc - ) V A((yt F )gi) nc . 

(67) 

It is clear from Eq. (|65|) that the non-constant term {{Y^Y) ! Ql) nc is a scalar. Further- 
more, from Eq. (I3TI) . we see that A is a scalar too. Since Yq has only scalar and vector 
components, the terms in Eq. ( JBTj) above have the same spinor structure. Thus, after 
taking the spinor trace, all these terms vanish. 

Putting together the above results, we have 

d 2 / 1 \ 1 r°° / 1 \ 2 / 

atrs ^_^_ 1 j r __/^„ (_) tr .^ (Y^Y Yi l) 

(68) 

One can perform the integration over a, since Yo((Y'Y)q1) c has no a dependence. 
We then have 

& tr '^ dJ ( y ^ 1 ) = ~2^ trsa - ( y o((^Sl)c) • (69) 

Since ((Y'Y)q1) c is a constant vector, when Iq acts from the left, the derivative terms 
in Yq vanish and, as a result, we have Yq((Y^Y)q1) c = (M — z6)((Yty)o'l) c . Since the 
first term in Eq. ( 1641) is a scalar, its contribution vanishes after taking the trace. We 
have then 

tr s0 ^yo(0^)oi) c 

= 2tr s <v (M - ifi) ^(^(0)^(0) - ^ q A /3 (0))(c q7q + irs a )( C/37/3 - irs^l 
= 2tr s «7^ 53 (yla (0)^(0) - 0ai4 y3 (O))(Mc a c j 37 Q! 7 y g + rs Q ,c ( gfj7 /3 - rs y3 c a fi7 a )l. 

(70) 
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To finish our calculation, we will use the fact that z and M are even functions of 
and thus any integral over k^ that involves an odd power of s M will vanish. For 
the spinor trace we use the relation tr s cr^jajp = ^i{^a^v(3 — b^pSua) and we finally 
obtain 



d 4 k d 2 / 1 " 

hm tr s o^v I Y , — ; — 1 



n (27r) 4 «-o h ^da 2 V WW J 

= "L(2^ ^ F ^ (0) - (71) 
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